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Abstract 

The finite size theory of met ast ability in a quartic potential is developed by the semiclassical 
path integral method. In the quantum regime, the relation between temperature and classical 
particle energy is found in terms of the first complete elliptic integral. At the sphaleron energy, 
the criterion which defines the extension of the quantum regime is recovered. Within the latter, 
the temperature effects on the fluctuation spectrum are evaluated by the functional determinants 
method and computed. The eigenvalue which causes metastability is determined as a function 
of size/temperature by solving a Lame equation. The ground state lifetime shows remarkable 
deviations with respect to the result of the infinite size theory. 

PACS numbers: 03.65.Sq - Semiclassical theories and applications. ll.10.Wx - Finite temperature field 
theory. 31. 15. Kb - Path integral methods. 74.50. +r - Tunneling phenomena 
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I. Introduction 



The theory of tunneling in a metastable potential has been a widely investigated topic 
after the seminal works by Langer and Coleman which provided the mathematical basis of 
the decay processes in statistical physics [l| and quantum field theory [2j. In the latter 
the decay rate of false ground states is suitably evaluated by the Euclidean path integral 
method in the semiclassical approximation 0, 4]: the classical particle paths are selected 
as the solutions of the Euler-Lagrange equations fulfilling some boundary conditions and, 
around such stationary points, the imaginary time action is expanded up to second order 
in the quantum fluctuations whose contribution is evaluated by means of the theory of the 
functional determinants. It is in the spectrum of the fluctuations that lie the origins of 
metastability jlj]. To be specific, take a particle of mass M moving in a nonlinear potential 
with frequency u and negative quartic parameter 5: 



Muj 2 2 5 4 

V(x) = — x - l x (!) 

Saying ±a are the positions of the potential maxima, S is given by 5 = Muj 2 /a 2 in units 
eVA~ A . Let's set a = 1A throughout the paper. The point x — is a classical ground 
state but quantum mechanically the particle can penetrate the hills and explore the abysses 
at \x\ > Vza. To see how it happens, let's focus on the positive a;— axis and consider the 
classical problem after performing a Wick rotation which maps the time from the real to 
the imaginary axis, t — > —ir, and turns the potential upside down as shown in Fig. [TJ 

Now the classical equation of motion admits a non trivial solution which, in the infinite 
size formalism, is known as the bounce. In the path integral language, the bounce is given by 
that path that starts at {x\ = 0, r = — oo), reaches the escape point at (x 2 = V2a, r = t ) 
and bounces back at (x\ = 0, r = +oo). The canonical bounce solution is therefore r— 
reversal invariant and consistent with the zero energy E motion (see Fig. [1]) between the 
turning points X\, Xi- In the latter the path velocity vanishes. As the bounce path has 
a maximum versus r the bounce path velocity has a node, hence the zero mode of the 
Schrodinger-like stability equation which governs the quantum fluctuations cannot be the 
ground state j|J . In other words, the bounce is a saddle and not a minimum for the Euclidean 
action. Thus, there must be a fluctuation whose eigenvalue is lower in energy than the zero 
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FIG. 1: (Color online) Reversed potential of Eq. JT]) (units meV). x\,X2 are the points in which 
the classical path velocity vanishes. Their location varies with E (in meV). to is in meV. 

mode eigenvalue. It is this negative eigenvalue which requires an analytic continuation in the 
Gaussian integral and ultimately leads to an imaginary square root fluctuation determinant. 
When multi bounce solutions are taken into account (due to multiple excursions between 
x\, X2) one finally gets an explicit formula for the semiclassical tunneling rate as given by 
the imaginary part of the ground state energy. The fact is that such tunneling rate does 
not depend on the temperature as the whole theory is based on the infinite size, or E = 0, 
formalism. The questions I want to deal with in this paper are the following: how is the 
quantum fluctuations spectrum affected by temperature effects in the quantum regime? 
And, accordingly, to which extent is the particle lifetime shortened? To get quantitative 
answers a finite size theory of metastability has to be developed. The focus is here on a non 



dissipative system 
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In Section II, I select the family of classical paths which makes stationary the Euclidean 
action thus generalizing the bounce solution for finite E. In Section III, the path integral 
of the metastable potential is presented together with the solution of the stability equation 
which governs the quantum fluctuations. The tunneling rate for the finite size theory is 
calculated in Section IV while the conclusions are drawn in Section V. 



II. Finite Size Classical Bounce 



In the imaginary time formalism 



9| the particle motion is classically allowed in — V(x) 



[Fig. CD) as the equation of motion reads: 
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Mx cl (r) = V'(x d {r)) 
where V means derivative with respect to x c i. Let's define 



(2) 



S 



and integrate Eq. (|2J). Then, one gets: 



r-r =±- [ dX 

(4) 

with integration constant E representing the classical energy associated to the particle 
motion. In Fig. [Tj, the potential barrier is sufficiently high to make the semiclassical approx- 
imation valid. The center of motion can be set at tq = with no loss of generality, then 
t G [-L/2, L/2] and L is the size of the system, that is the period for a particle excursion to 
the abyss and back to the starting position. Then the size represents the finite time required 

c 4 

for this journey to occur. For negative energies, — ^f- < E < 0, there are two turning points 
x i (Xi) an d x 2 {xi) (0 < < x 2 ) which define the bounds for the particle excursion. They 
are given by: 



xi =y — 2 — ;x 2 =y — g — (5) 

Thus the length of the generalized bounce, X2 — Xi> depends on S and shrinks from 1 
(i? = 0) to {E = —V(x = a)) [h]]. This marks an essential difference with respect to the 
bistable 4 model ll| in which the (anti)instantons interpolates between the two potential 



minima thereby covering a distance that does not depend on E. To integrate Eq. (0J I use 
the result 
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X2 



dx 



Xcl 



M {x 2 -xi){xl-x 2 ) 



—F((,m) 

X2 



( = arcsin 

2 



'Xl -Xciirf 



xl~xl 



i - xilxi 



(6) 



with F((,m) being the elliptic integral of the first kind with amplitude ( and modulus 
m. Then, after setting Xciijo) = X2, from Eqs. (J3j) , (J4j) , ([6]) I get the generalized bounce for 
the finite size theory: 



x c i{t) = a\ / 2x2dn(w,m) 

ZU = X2^(t - T ) 



(7) 



where dn(zu,m) is the Jacobi delta-amplitude defined in Eq. (IA.1I) . In the E = 
limit, m 2 = 1, hence one recovers the bounce solution of the infinite size theory, x c i(r) = 
a\f2 cosh -1 (ro). Fig. [2] plots some paths of the family in Eq. (J7J): the shape of the paths 
is not essentially modified (with respect to the E ~ path) up to energies of order 
\E\ ~ 0.1V(x = a). Above this value the bounce progressively shrinks and becomes a 

131 ] defined by E sph = \E\ — V(x = a). 



point-like object at the sphaleron energy 



dn(-cu,m) is even function of w in the period [— K(m), K(m)] determined by the first 
complete elliptic integral K(m) = F(7r/2,m) [14]. Then, imposing the boundary conditions 
on the t— range, I get from Eq. (CO) (with r = 0): 



00L 



2K(m) 

X2 



(8) 



Eq. ([H]) expresses the fundamental relation between the size L and the particle energy E 
which is hidden in the modulus according to Eqs. ([3]), © and the last of Eq. Mapping L 
onto the temperature axis, L = h/ (KbT*) in the spirit of the thermodynamic Bethe Ansatz 
[15I ]. one obtains the link between E and the temperature T* at which the particle motion 
takes place: 



K B T* 



huJ X2 
2 K(m) 



(9) 
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For E = — V(x = a), Kim = 0) = tt/2 hence, Eq. (Q leads to the Goldanskii criterion 
[3] for the transition between quantum and activated regime: KbT^^ = Huj / (y/2~7r) . T^ ax 
is thus the maximum temperature at which the tunneling occurs. Computation of Eq. (IS]) 
shows that L(|J5|) is a monotonically decreasing function below the sphaleron energy. This 
confirms, on general grounds that the transition at T^ ax is expected to be a smooth 
crossover as proposed long time ago jf], [ijj]. The physical origin of the smooth change lies 
in the fact that the bounce is flexible and continuously adapts its shape to the periodic and 



changeable (with E) boundary conditions. Thus the physical 



that encountered in studying the bistable 4 potential [H 
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picture differs very much from 



20] where the (anti)instantons 



have to fulfill antiperiodic boundary conditions which are simply imposed by the potential 
structure. As a consequence the quantum/activated crossover for the bistable 4 potential 



is in fact a sharp transition 
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The same conclusion regarding the character of the transition at T^ ax for the metastable 
<fi potential may be drawn by a direct evaluation of the classical action A[x c ;] which, for the 
finite size bounce, reads 



A[x d ] = 2V2M j dx^E + V(x) -E-L 



(10) 



Using Eqs. ©, (CJ), ([H]) one can monitor the smooth E (or T*) dependence of A[x c j\ and 



its derivative 
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For E — » 0, Eq. ffTUj) leads to the infinite size result 



A[x d }/h -> 4M 2 w 3 /(3h8) (11) 

The fact that A[x c j] oc 1/8 represents the motivation for the semiclassical approach to 
the quantum tunneling. 

III. Semiclassical Path Integral 

The particle path can be written as a sum of the classical background and the quantum 
fluctuation, x(r) = x c i{t) + 77 (r). Then, the semiclassical space-time particle propagator 
between the positions Xi and Xf in the imaginary time L is given, in quadratic approximation, 
by 
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FIG. 2: (Color online) Shape of the bounce solution for various classical energies E. L is in sec 
and T* in K. 



< Xf\xi >l= exp 

-L/2 



-^A[x d } 



Dr] exp 



dr 



M ( 



L/2 



\nr) + -V"(x cl (r)W(r) 



M 



V"{x d {r)) = Mu 2 (l - 6xldn 2 (w,m)) 



(12) 



As Xi and Xf coincide in our model, Eq. fll2jl is the single bounce contribution Z\ to the 
total partition function which determines the tunneling rate. Let's evaluate it. The measure 
Drj of the path integration is given through the coefficients q n of the fluctuation expansion 
in a series of ortonormal components rj n (r): 



/ °° PCX. 

oo 
n=0 



fan 



2nh/M 



(13) 



The normalization constant K accounts for the Jacobian in the transformation to the 
normal mode expansion and the components r] n (r) [22J are the eigenstates (whose eigenvalues 
are denoted by e n ) of the Schrodinger-like equation: 



d 2 V"{x d {r)) 



r} n {r) = e n 7] n (T) 



(14) 



dr 2 ' M 

Supposed to have solved Eq. (fl4"|) (see below), after performing Gaussian integrations 
over the directions one formally obtains the quantum fluctuations factor: 
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Drjexp 



h 



X(Det[0] 



-1/2 



Det[6] = 11 o 

n=0 



E n - 0= -dt + 



2 , V"(X CI ) 



M 



(15) 



There is however a trouble in Eq. ({15]) due to the zero eigenvalue £o which breaks the 
Gaussian approximation making Det[0\" 1 / 2 divergent. It can be easily seen [f]] that the 
zero mode is 23( the path velocity x c i(r) which solves the Euler-Lagrange problem and the 
homogeneous differential equation associated to the Eq. ( fl4"l) . The physical origin of the 
zero mode lies in the fact that the center of motion r of the periodic bounce solutions can 
be placed everywhere in the range [-L/2, L/2}. This mirrors the r-translational invariance 
of the system. Then the zero eigenvalue can be extracted from Det[0] and, resorting to 
the integration over oi in Eq. (fl3l) . it can be replaced in Eq. ( fl"5l) according to the recipe: 

1-1/2 _ 



A[x cl \ 
2nh - 



24] 



A. Functional Determinant 
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27 



28, 29, 3 



Now I proceed invoking the theory of the functional determinants 
to compute the whole fluctuation contribution embodied in the regularized determinant 
Det R [6] such that Det[6) = e ■ Det R [0]. The form of Det R [6] depends on the type 
of boundary conditions fulfilled by the quantum fluctuations rj n {j) and their derivatives. 
To establish it, remember that x c i(t) is a quantum fluctuation whose periodicity can be 
promptly checked by deriving Eq. ([7|): 



%ci( T ) = — aoj\/2rri 2 'x^sn{w, m)cn(w, m) 

Xd( T ) = —cllo 2 \f2m 2 y^dn{vj ) m) [l — 2sn 2 [w, m)] 

(16) 

with the Jacobian elliptic functions sn(zu,m) and cn(zu,m) defined in Eq. ( 1A.1I) . Then, 
for any two points such that w 2 = vj\^ 2K(m), one infers that: x c i{vj2) = x c i{vji) and 
Xd{^2) = x c i{w\). Thus, periodic boundary conditions (PBC) apply to our problem. 
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In fact, being a product over an infinite number of eigenvalues whose modulus is larger 
than one, Det R [0] is (exponentially) divergent in the L — > oo limit whereas ratios of func- 



tional determinants are finite and physically meaningful both in value and sign [25|, |31 |. 
Therefore Det R [0] has to be normalized over the harmonic oscillator determinant Det[h] 
with h = —d 2 + uj 2 . For any two points separated (as above) by 2K(m), the latter corre- 
sponding to the period L along the r— axis, the determinants ratio is formed in the case of 



PBC by 



33|: 



Det[h) = -4sinh 2 (cjL/2) 
Det R [0] fx{w 2 )-h{w x ) 



</o|/o> /o(^i)W(/o,/i) 
, dx d (r) 

J 1 ~ — ^ 

urn 

h = x d (r) (17) 

where f , fi are independent solutions of the homogeneous equation associated to 
Eq. ffl^l) . W(fo, fi) is their Wronskian and < fo\fo > is the squared norm which can be 
computed by the first in Eq. f|T6|) . 

Working out the analytical calculation for the second in Eq. (jTTl) . I find: 



Det R [6] (2 - m 2 ) 7 ' 2 



< k\k > a 2 u 3 m 4 



— ; - 2 K { m ) 

1 — m z 



(18) 



where E(7i/2,m) is the complete elliptic integral of the second kind. As the squared 
norm has dimension [ooa 2 ], Det R [0] correctly carries the dimension [u)~ 2 ] consistently with 
the fact that Det[0] is dimensionless. 



In the E — > (L — > oo) limit, m 2 — > 1 and K{m) ~ ln(4/\/l — m 2 ). Moreover, < 



/olio >~ f^a 2 - 



Then, from Eqs. (|T7j) . ( TTBl . in the same limit, I get 



Det[h] _ , 



12oj 2 (19) 



Det R [0] 

Thus recovering the well known result of the infinite size bounce theory. In fact this 
occurs quite away from the zero limit as clearly shown in Fig. |3] where the determinants 



9 



Q 4 

o 

92 




- 4 



12 



FIG. 3: (Color online) Normalized determinants ratio versus T*(K). The inset plots the regularized 
fluctuation determinant (units meV~ 2 ) of Eq. (fT8l) . 

ratio (normalized over — 12uo 2 ) is plotted against T* obtained from Eq. ([9]). As u = 5meV, 
T^ ax is set here at ~ 13K although it may be much larger in real systems 32j. The result 
of Eq. (fT9|) is already achieved at T* ~ 7K, that is about the value at which the finite size 
bounce fits the infinite size bounce (see Fig. [2]). This looks remarkable as it suggests that 
quantum fluctuations and classical path start to experience the finite size effects starting 
from the same temperature. In the T* — > limit, computation of Eq. (1TB]) becomes time 
consuming in order to control the elliptic integrals and reproduce the exponential divergence 
due to the (1 — m 2 ) -1 term. 

The negative sign in Eq. (fT9l) is crucial and occurs for any T*. Then, the square root of 
the determinant ratio (which appears in the partition function) is imaginary. As stated at 
the beginning and made clear in Fig. [2] the bounce has always a maximum hence, the zero 
mode is not the ground state: there is a negative eigenvalue £_i 34] embedded in Det R [0] 
which causes a finite lifetime for the particle placed at x = x\. While such eigenvalue is 
easily found = — 3u 2 ) in the infinite size theory as the operator in the l.h.s. of Eq. (fT4"j) 
becomes of Rosen- Morse type 35j , here I am facing a non trivial problem: that to determine 
£_i at finite E and establish its effect in the temperature dependence of the particle lifetime. 



B. Lame equation 

To pursue this goal, insert the last of Eqs. (|T2|) in the stability equation (fT4|) which easily 
transforms into: 
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d 2 

Vn(r) = [1(1 + l)m 2 sn 2 (zu,m) + A n ]rin(r) 



dw 2 

A n 9 1(1 ~\~ i) ~, 

xi (x2uj) 2 

1(1 + 1) = 6 (20) 

This is the Lame equation in the Jacobian form [3^] that, for given / and m, admits 
periodic solutions (which can be expanded in infinite series) for an infinite sequence of char- 
acteristic A n values. The continuum of the fluctuation spectrum stems from this sequence. 
However, for positive and integer /, (in the case at hand, 1=2), the first 21 + 1 solutions 
of Eq. fl20|) are not infinite series but polynomials in the Jacobi elliptic functions with real 
period 2K(m) or AK(m) [3?]]. 2K(m) plays the role of the lattice constant being the period 
of the potential in the Schrodinger like stability equation. It is among these five polynomial 
solutions corresponding to the eigenvalues A n that one has to pick up the unstable fluctua- 
tion eigenstate. Two solutions out of five have to be discarded because they do not fulfill the 
periodicity conditions required for the fluctuation components: rj n (w) = rj n (w =F 2K(m)). 
The three (unnormalized) good solutions are 



77_i oc sn 2 (w,m) -==== 

1 + m 2 — v m 4 — m 2 + 1 

?7o oc sn(zu,m)cn(zu,m) 

rji oc sn 2 (w,m) -=== 

1 + m? + v m 4 — m? + 1 

(21) 

and, using Eqs. fl20l . I get the corresponding energy bands e n as functions of the modulus 

m: 



2 



S-i = —u 2 I 1 H -Vm 4 — m 2 + 1 

2 — m 2 



e = 



2 



e l= -u 2 1-- Vm 4 - m 2 + 1 (22) 

\ 2 — m z I 

The zero mode is thus consistently recovered also for the finite size theory with r] oc x c \ 



23j. r/i lies in the continuum with ei — > uj 2 in the infinite size limit. e_! is the object of our 
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FIG. 4: (Color online) Ratio (in meV~ A ) between the regularized fluctuation determinant and the 
ground state eigenvalue versus temperature. The inset plots separately the T* dependence of the 
g.s. eigenvalue (in meV 2 ). 
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FIG. 5: (Color online) Tunneling rate (in meV) computed from Eq. f)23|) versus temperature. The 
rate of the infinite size theory is plotted for comparison. The inset shows the lifetime (in sec) 
against the size L. 

focus. Through Eqs. ([5]), ([6]) and ([9]), the T* dependence of e_i is computed as shown in 
the inset of Fig. HJ the value —3a; 2 of the infinite size theory is mantained up to T* ~ 7K 
while softens significantly at larger T*. This strongly affects the overall behavior of 
Det R [0] as I point out by plotting the ratio Det R [0}/e-i on a linear scale. 

IV. Decay Rate 

All the ingredients are now available to compute Eq. fTl2l) and derive the total partition 
function Zt- Note that also the solution of Eq. ([2]), x c \ = contributes with a term Zh to 
Z T . Being Z h = H\Det[h] | -1 ^ 2 , we now fully understand on physical grounds the need to 
normalize Det R [0] over Det[h] as described above. In fact, one has to account for all multiple 
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excursions to and from the abyss which is equivalent to sum over an infinite number of (non 
interacting) bounce contributions as Z\ in Eq. (j!2p . The final result is: Zt = Zh exp(Zi/Zh)- 
The tunneling rate T(L) is thus obtained from the imaginary part of the exponent in Zt 
(Zi is purely imaginary) through the Feynman-Kac formula [5|. The general expression for 
a finite size system is: 



(23) 

with Det RR [0] = Det R [0]/e-i being negative as shown in Fig. HI Eq. (123]) is plotted in 
Fig. [5] against T* together with the constant tunneling rate of the infinite size theory. T(T*) 
grows fast as the length of the bounce shrinks and the quantum fluctuations spectrum 
softens. Close to T^ ax , the tunneling rate has increased by an order of magnitude with 
respect to the infinite size result but the decay width remains smaller than the fundamental 
oscillator energy being T(T* = 12.98K)/uj ~ 6 • 10~ 3 . The inset in Fig. [5] plots the lifetime 
of the false vacuum versus the system size. 



r(L) = n 



A[x d ] +E-L 



Det[h] 



2nh \ Det RR [6)\e^ 



exp 



V. Conclusion 



This work presents a study of the finite size effects on a metastable quartic nonlinear 
potential and, as such, is complementary to a previous investigation The semiclassical 
path integral method, implemented by the theory of the functional determinants, builds 
the framework to analyse the tunneling in the metastable system. The description relies 
on the properties of the elliptic functions which permit to monitor the evolution of the 
classical bounce versus the system size. After defining the temperature range within which 
the tunneling occurs, I have computed the temperature effects on the quantum fluctuation 
spectrum. In particular, the negative eigenvalue which causes metastability has been studied 
in great detail solving a Lame type equation. These results, new in the literature, permit to 
compute the physical properties of the false ground state. Specifically I have shown that its 
lifetime has a remarkable size/temperature dependence inside the quantum regime. While 
the latter can persist up to T^ ax of order 100K in real systems, the presented quantitative 
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method may also be of practical interest provided the potential parameters are adapted to 
specific cases. 

APPENDIX 

The Jacobian elliptic functions used throughout the paper are related to the amplitude 
( in Eq. (jH]) by the definitions 



sn(w, to) = sin£ 
cn(zcr, m) = cos( 

dn(w,m) = \J 1 — m 2 sin 2 ( (A.l) 

In the computation of the bounce and its time derivative, the following representations 
as trigonometric series have been used [l2| 

2tt ^ q n - x l 2 . /(2n-lW\ 

sn{w,m) = ; — - > — -sin — - — 

V ' ; mK(m) ^ 1 - q 2n ~ x V 2K(m) ) 

2tt q n - l l 2 f (2n- l)irw \ 

mK(m) 1 + q 2n ~ l C \ 2K(m) ) 

ix 2ix q n ( nnzu 



. . 71 Z.7I X q / 

dn(w, m) = ——. — - + — - — - > — cos -—, — r 

v ' ' 2K(m) K(m) ^ 1 + q 2n \K(m). 

q = exp(—7rK'(rn)/K(m)) 
K'{m) = F(tt/2,TO) 

m 2 = 1-m 2 (A.2) 

with w defined in Eq. (J7J). Numerical convergence is achieved by taking a cutoff n max ~ 40 
in the Fourier series. 
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related to the instanton and bounce trajectories have zero and one node respectively. 
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